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Abstract 

We consider singular perturbation elliptic problems depending on a parameter e such that, for 
e = the boundary conditions are not adapted to the equation (they do not satisfy the Shapiro 
- Lopatinskii condition). The limit only holds in very abstract spaces out of distribution theory 
involving complexification and non-local phenomena. We give a very elementary model problem 
showing the main features of the limit process, as well as a heuristic integral procedure for obtain- 
ing a description of the solutions for small e. Such kind of problems appear in thin shell theory 
when the middle surface is elliptic and the shell is fixed by a part of the boundary and free by the 
rest. 

1 Introduction 

The main purpose of this paper is to give general ideas on a kind of singular perturbations arising 
in thin shell theory when the middle surface is elliptic and the shell is fixed by a part of the boundary 
and free by the rest as well as an integral heuristic procedure reducing them to simpler problems. The 
system depends drastically on the parameter e equal to the relative thickness of the shell. It appears 
that the "limit problem" for e = is highly ill-posed. Indeed, the boundary conditions on the free 
boundary are not "adapted" to the system of equations; they do not satisfy the Shapiro - Lopatinskii 
(SL hereafter) condition. Roughly speaking, this amounts to some kind of "transparency" of the 
boundary conditions, which allow some kind of locally indeterminate oscillations along the boundary, 
exponentially decreasing inside the domain. This pathological behavior is only concerned with e = 0. 
In fact, for e > the problem is "classical". When e is positive but small, the "determinacy" of the 
oscillations only holds with the help of boundary conditions on other boundaries, as well as the small 
terms coming from e > 0. 

In such kind of situations, the limit problem has no solution within classical theory of partial 
differential equations, which is uses distribution theory. It is sometimes possible to prove the con- 
vergence of the solutions u £ towards some limit u°, but this "limit solution" and the topology of the 
convergence are concerned with abstract spaces not included in the distribution space. 

After recalling the SL condition (section [2]), we give in section [3] a very simple example of such 
a perturbation problem. The geometty of the domain (an infinite strip) allows explicit treatment by 
Fourier transform in the longitudinal direction. The inverse Fourier transform within distribution the- 
ory is only possible for e > 0, whereas for e = it is only possible in the framework of analytic 
functionals (higly singular and not enjoying localization properties). This example shows the promi- 
nent role of components with high frequency; for small e, the "smooth parts" (i. e. with small |£|) of 
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the solutions may be neglected with respect to "singular ones" (i. e. with large |£|). We also recall 
an example of elliptic Cauchy problem (in fact Hadamard's counter-example) which exhibits some 
relation with the limit problem. 

In sectionlU we report the heuristic procedure of [EgMeSa07 ]. In this latter article, we addressed a 



more complicated problem including a variational structure, somewhat analogous to the shell problem, 
but simpler, as concerning an equation instead of a system. It is shown that the limit problem contains 
in particular an elliptic Cauchy problem. This problem was handled in both a rigorous (very abstract) 
framework and using a heuristic procedure for exhibiting the structure of the solutions with very small 
e. The reasons why the solution goes out of the distibution space as e goes to are then evident. In 
section [4] we present a simplified version of the heuristic procedure involving only the essential facts 
of the approximation, which are very much analogous to the method of construction of a parametrix 
in elliptic problems llTa8TTl . (EgSc97j: 

-Only principal (with higher differentiation order) terms are taken into account. 

-Locally, the coefficients are considered to be constant, their values being frozen at the corre- 
sponding points. 

-After Fourier transform (x — > £)> terms with small £ are neglected with respect to those with 
larger £ (which amounts to taking into account singular parts of the solutions while neglecting smoother 
ones). We note that this approximation, aside with the two previous ones, lead to some kind of "local 
Fourier transform" which we shall use freely in the sequel. 

Another important ingredient of the heuristics is a previous drastic restriction of the space where 
the variational problem is handled. In order to search for the minimum of energy, we only take 
into account functions such that the energy of the limit problem is very small. This is done using 
a boundary layer method within the previous approximations, i. e. for large |£|. This leads to an 
approximate simpler formulation of the problem for small e, where it is apparent that the limit problem 
involves a smoothing operator and cannot have a solution within distribution theory. 

It should prove useful to give an example of a sequence of functions converging to an analytical 
functional (but going out of the distribution space, then leading to a "complexification" phenomenon). 
It is known ([Sc50], [GeCh64]) that (direct and inverse) Fourier transform within distribution theory 
is only possible for temperate distributions, not allowing functions with exponential growth at infinity. 
The space of (direct or inverse) Fourier transform of general distributions is noted Z' . It is a space 
of analytical functionals: the corresponding test functions are analytical rapidly decreasing functions, 
forming the space Z. 

Let us consider the (non temperate) distribution (or function) u{£) = cosh(^). The sequence 



u A (0 



cosh(0 if |£| < A, 
otherwise 



converges to u in the distribution sense as A goes to infinity. The inverse Fourier transforms u x (x) 
converge in Z' to the analytical functionnal u(x). The functions u A (£) are tempered and their inverse 
Fourier transforms are easily computed by hand. It appeals that for large A 



e A 1 



u (x) ~ ^-(cosfAx) + xsin(Ax)). 

2tt 1 + x z 



It is then apparent that u (x) consists of a "nearly periodic" function with period tending to zero 
along with 1/A, multiplied by an "envelop" defined by jt^i and by the factor Moreover, it should 
be noted that the amplitude is exponentially large with respect to the inverse of the period. It is 
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Figure 1.1: Normal displacement for e = 10 3 on the left and for e = 10 5 on the right 

then apparent that the limit is an "extremely singular" function as the "graph" fills the entire plane. 
Moreover, it is clear (and may be rigorously proved [EgMeSa07 ] that the sequence u x goes out of the 



distribution space everywhere, not only in the vicinity of x = as is suggested by the formal inverse 
Fourier transform of cosh(£) = ^^=0 (In)!' which is 



U{X) = ^=°j2^ 62n{x) > 

apparently a singularity "of order infinity" at the origin. This fact constitutes an example of the prop- 
erty that elements of Z' can only be tested with analytical test functions, then not enjoying localization 
properties. 

The motivation for studying that kind of problems comes from shell theory, see |SaHu Sa97l , 
[BeMiSa08]. It appears that when the middle surface is elliptic (both principal curvatures have same 
sign) and is fixed by a part Tq of the boundary and free by the rest Ti, the "limit problem" as the 
thickness e tends to zero is elliptic, with boundary conditions satisfying SL on To, and boundary not 
satisfying SL on Fi. Without going into details, which may be found in [MeSa06], [MeSaHuSa07], 



|EgMeSa07J and |EgMeSa09], we show numerical computations taken from [BeMiSa08 ] of the nor- 
mal displacement for e = 10 -3 and e = 10 -5 (figures [TTTI on the left and on the right respectively) 
when the shell is acted upon by a normal density of forces on a rectangular region of the plane of 
parameters. The most important feature is constituted by large oscillations nearby the free boundary 
Y\. It is apparent that, when passing from e = 10~ 3 to e = 10~ 5 , the amplitude of the oscillations 
grows from 0.001 to 0.01. The singularities produced by the jump of the applied forces inside the 
domain is still apparent for e = 10~ 3 , not for e = 10 -5 , where only oscillations along the boundary 
are visible. Moreover, the number of such oscillations pass from nearly 3 for e = 10~ 3 to nearly 5 for 
e = 10~ 5 and is then nearly proportional to log(l/e). We shall see that all these features agree with 
our theory. 

2 The Shapiro - Lopatinskii condition for boundary conditions of el- 
liptic equations 



In this section, we recall some properties of elliptic PDE, see [ AgDoNi59 ] and [EgSc97 ] for more 
details. 
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We consider a PDE of the form 

P(x,d a )u = f{x) (2.1) 

Where x = [x\, X2) and d a = d/dx a , a = 1,2, and P is a polynomial of degree 2m in d a . Let Po 
be the "principal part", i. e., the terms of higher order. The equation is said to be elliptic at x if the 
homogeneous polynomial of degree 2m in £ a : 

P o (x,-i^) = (2.2) 

has no solution £ = (£1, £2) 7^ (0, 0) with real £ a . When the coefficients are real (this is the only case 
that we shall consider) this implies that the degree is even (this is the reason why we denoted it by 
2m). The left hand side of (12.21 ) is said to be the "principal symbol"; the "symbol" is obtained in an 
analogous way taking the whole P instead of the principal part Pq. We note that replacing d/dx a by 
—i^a in Pq amounts to taking formally the Fourier transform x — > £ for the homogeneous equation 
with constant coefficients obtained by discarding the lower order terms and freezing the coefficients 
at x. Obviously, ellipticity on a domain Q, is defined as elliptic at any x £ Q. 

It is worthwhile mentioning that ellipticity amounts to non - existence of "travelling waves" of the 
form 

e~^ x (2.3) 

for the equation obtained after discarding lower order terms and freezing coefficients. Here "travel- 
ling" amounts to "with real note that solutions as (12.31 ) with non real £ are necessarily exponentially 
growing or decaying (in modulus) in some direction. Moreover, when a solution of the form (12.31 ) ex- 
ists (with £ either real or not), it also exists for c£ with any c. In a heuristic framework, we may suppose 
that |£| is very large; this justify to discard lower order terms (= of lower degree in |£|). In the same 
(heuristic) order of ideas, freezing the coefficients allows to consider "local solutions". This amounts 
to multiply the solutions by a "cutoff" function 9{x) or equivalently taking the convolution of the 
Fourier transform with #(£), which do not modify the behavior for large £. Microlocal analysis gives 
a rigorous sense to that heuristics. It then appears that local singularities of a solution u (associated 
with behavior of the Fourier transform for large |£|) cannot occur in elliptic equations unless they are 
controlled by the (Fourier transform of the) right hand side /. This gives a "heuristic proof" of the 
classical property that local solutions of elliptic equations are rigorously associated with singularities 
of/. 

What happens with solutions near the boundary? Local Fourier transform is no longer possible, 
but, after rectification of the boundary in the neighborhood of a point, we may perform a tangential 
Fourier transform. If, for instance, the considered part of the boundary is on the axis x\ and the domain 
is on the side X2 > 0, taking only higher order terms and frozen coefficients, we have solutions of the 
form (12.31 ) with real £1 (as coming from the Fourier transform) and non - real £2. The dependence in 
X2 is immediately obtained by solving an ODE with constant coefficients. Obviously, the solutions 
are exponentially growing or decreasing for x-i > 0. As the coefficients are real, there are precisely 
m (linearly independent) growing and m decreasing (in the case of multiple roots, dependence in X2 
of the form X2e K ' X2 and analogous also occur). Roughly speaking, there are solutions of the form: 

^2c k e-^ lXl e XkX2 (2.4) 

k 

with real £1 and Re(X) 7^ (here k is running from 1 to 2m) . Boundary conditions on X2 = 
should control solutions with Re(X) < 0, i. e., exponentially decreasing inside the domain, whereas 
exponentially growing ones should be controlled "by the equation in the rest of the domain and the 
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boundary conditions on the other parts of the boundary". In other words, "good boundary conditions" 
should determine, (within our approximation of the half plane and frozen coefficients) the solutions of 
the equation of the form (12.41 ) with Re(X) < 0. Obviously, the number of such boundary conditions 
is to. A set of m boundary conditions enjoying the above property is said to satisfy the Shapiro - 
Lopatinskii condition. There are several equivalent specific definitions of it. We shall mainly use the 
following one: 

Definition 2.1. Let P be elliptic at a point O of the boundary. A set of m boundary conditions 
Bj(x, d a ) = gj(x), j = 1, ...to is said to satisfy the SL condition at O when, after a local change to 
new coordinates with origin at O and axis x\ tangent to the boundary, taking only the higher order 
terms and coefficients frozen at O in the equation and the boundary conditions, the solutions of the 
form ri2.4D with Re(X) < obtained by formal tangential Fourier transform are well defined by the 
boundary conditions. 

Remark 2.1. The above definition should be understood in the sense of formal solution for any given 
(real and non-zero) £i. The SL condition is not concerned with solutions in certain spaces. It is purely 
algebraic, and concerns to conditions imposed to the m (decreasing with x 2 ) linearly independent 
solutions of the ODE obtained from Pq by formal tangential Fourier transform. This also amounts to 
saying that imposing the boundary conditions equal to zero, the considered solutions must vanish. In 
fact, the SL condition amounts to non-vanishing of a certain determinant, and as so it is generically 
satisfied: conditions not satisfying it are rarely encountered. In particular, in "well-behaved prob- 
lems", when coerciveness on appropriate spaces is proved, the SL condition is not usually checked. 
It should also be noted that the SL condition is independent of a change of variables, and, in most 
cases, the change is trivial. On the other hand, there are also definitions of the SL condition without 
change of variables. Last, it should also be noted that the SL condition has nothing to do with lower 
order terms and the right hand side of the boundary conditions ( as ellipticity is only concerned with 
the principal symbol); it is merely a condition of adequation of the principal part of the boundary 
operators to the principal part of the equation. 

Let us consider, as an exercise, examples for the laplacian. 

p = -d\ - d 2 2 (2.5) 

The principal symbol is £f + £| so tnat the equation is elliptic of order 2, then to = 1. "Good 
boundary conditions" are in number of 1. 

Let us try the boundary condition (Dirichlet): 

u = 0. (2.6) 

Taking any point of the boundary and {x\,xz) with origin at that point, tangent and normal to the 
boundary respectively, the equation is the same as in the initial variables, and formal tangential Fourier 
transform gives 

(£i 2 -d 2 2 )n(£i,x 2 ) = (2.7) 

and the solutions are 

u{ti,x 2 ) = d(6)el^' 2 + C 2 (6)e lei|a:2 . (2.8) 
Taking only the exponentially decreasing for x 2 > we only have 

u(&,x 2 ) = Ci^e-I^ 2 . (2.9) 
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Now, imposing the "tangential Fourier transform" of (12.61 ): 



6(6,0) =0, (2.10) 

we see that it vanish identically. Then, the Dirichlet boundary condition satisfies the SL condition for 
the laplacian. 

The case of the Neumann boundary condition for the laplacian 

^ = 0. (2.11) 

on 

is analogous. (Note also that the Fourier condition (|^) + au = g is the same, as only the higher 
order terms are taken in consideration). Proceeding as before, we have, instead of (12.101 ): 

02^1,0) = -|6|Ci = 0, (2.12) 



which also gives Ci(£i) = and then u = 0. Thus, (12.101 ) satisfies SL for (I2.5I ). 
Oppositely, the boundary condition: 

(d s -id n )u = 0, (2.13) 

where s and n denote the arc of the boundary and the normal, does not satisfy the SL condition for the 
laplacian. Indeed, taking the new local axes, s and n become x\ and X2, and after tangential Fourier 
transform: 

Kl-»W(ft,0) = 0, (2.14) 

which, applied to (12.91 ) becomes: 

H6 + i|6l)Ci(6) = o. (2.15) 

we then see that C\ (£i) vanishes for negative £i, but is arbitrary for positive £j. In fact , the boundary 
condition (12.131 ) is "transparent" for solutions of the form ( 12.91 ) with positive £i. 

Remark 2.2. As it is apparent in the last example, when the SL condition is not satisfied, there is 
some kind of "local non-uniqueness", where "local" recalls that only higher order terms are taken in 
consideration, and the coefficients are frozen at the considered point of the boundary. 

The SL condition appears as some previous condition for solving elliptic problems. It is apparent 
that some pathology is involved at points of the boundary where it is not satisfied. 

Let us mention, before closing this section, that the boundary conditions may be different on 
different parts of the boundary specially on different connected components of it (when there are 
points of junction of the various regions, usually singularities appear at that points). 



3 An explicit perturbation problem where the SL condition is not satis- 
fied on a part of the boundary of the limit problem 

Let SI be the strip (— oo, +oo) x (0, 1) of the (x, y) plan. We denote by To and T\ the boundaries 
y = and y = 1 respectively. We then consider the boundary value problem depending on the 
parameter e: 

Au £ = on Q 

n e = 0onr (3.16) 

d x u + (i + e 2 )d y u = ip on Ti 
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where tp is the data of the problem. It is a given function of x, that we shall suppose sufficiently 
smooth, tending to at infinity. We shall solve it by x — > £ Fourier transform; it is easily seen that we 
also have automatically u — > for x — > oo, which may be added to the boundary conditions. 

The boundary condition on Tq is the Dirichlet one, which satisfies SL for the laplacian. Oppositely, 
the boundary condition on Ti satisfies it for e > (this is easily checked), not at the limit e = (see 
the end of the previous section). The problem is to solve for e > and to study the behavior for e 
going to zero. 

Denoting by "the x — > £ Fourier transform, u e is defined on the same $7 domain, but of the (£, y) 
plane. The solutions of the (transform of) equation and the boundary condition on To are of the form 

u £ {^y) =a(£)sinh(£y) (3.17) 

where a denotes an unknown function to be determined with the boundary condition on Ti. It will 
prove useful to write the solution under the form 

smn(4) 

for the new unknown /3 £ (£), which is the transform of the trace u £ (x, 0). Imposing the Fourier trans- 
form of the boundary condition on F± we have: 

- + + e 2 )^$-P £ m = 0(0- (3-19) 
smn(4J 

So that: 

= 7 ■ (3-20) 

—i£ ( 1 — coth(f ) ) + e 2 i coth(£) 



In order to study this function, we should keep in mind that the expression (1 — coth(£)) decays for 
£ — ► +cxd as 2e -2 ^. Then, at the limit e = we have 



= wr^m - (3 - 21) 

- coth(O) 



For £ —>■ +oo this function behaves as 



00^ ra 2 M e 2«. ( 3.22) 



This shows (unless in the case of very special data (p with very fast decaying Fourier transform) that 
is not a tempered distribution, and the inverse Fourier transform is an analytical function in Z'. 
Nevertheless, for e > 0, is "well-behaved" for £ — ► +oo as 



P £ (0 ~ (3-23) 



This specific behavior depends on that of |-, so that in most cases will be decreasing, but multiplied 

by the factor e~ 2 . When e > (small but not 0) is fixed, /3 e (£) is approximatively given by (13.211 ) 
for "finite" £ and by (13.231 ) for £ going to +oo. It is easily seen that the sup in modulus of |/3 £ (£)| is 
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located in the region where both terms in the denominator of the right hand side of (13.201 ) are of the 
same order (so that no one of them may ble neglected). This gives 

£ = 0{log(l/e)). (3.24) 

It appears that /3 £ (£) consists mainly of Fourier components which tend to infinity algebraically 
as e goes to zero with £ tending to infinite "slowly" as in (13.241) . this is somewhat analogous to the 
example, given in the introduction, of a sequence of functions converging to an analytical functional. 

Coming back to (13.181 ). the main properties of the behavior of u e (x, 1) may be thrown: 

- The trace u e (x, 1) = f3 e {x) on the boundary Y\ which bears the "pathological boundary condi- 
tion" mainly consists of large oscillations with wave length l/log(l/e) (which tends to very slowly 
as e — > 0). The amplitude of that oscillations grows nearly as e~ 2 . The limit e — > does not exist in 
distribution theory; it constitutes a complexification process. 

- Out of the trace on Y\, (i. e. for < y < 1), the behavior is analogous, but of lower amplitude, 
which is exponentially decreasing going away of Y\. We recover properties of the non-uniqueness 
associated with the failed SL condition. 

Before concluding this section, we would like to show some analogy between the previous limit 
problem and the Cauchy elliptic problem, which is a classical example of ill-posed problem, without 
solution in general. 

We consider the same domain £1 as before, but we now impose two boundary conditions on Tq 
and no condition on T\. Namely 

Av = on n 

v = ip on T (3.25) 
d y v = on To 

Taking as above the x — > £ Fourier transform, it follows immediately that 

£(£,y)=*0(£)eosh(£y). (3.26) 

Where it is apparent that the behavior for £ — > oo is exponentially growing (unless in the case when 
decays faster than e - '^) so that it is not tempered and the inverse Fourier transform does not 
exist within distribution theory. 



4 A model variational sensitive singular perturbation, [EgMeSa07] 

4.1 Setting of the problem 

Let f2 be a two dimensional compact manifold with smooth (of C°° class) boundary dQ = TqUFi 
of the variable x = (xi, 22), where Tq and Ti are disjoint; they are one - dimensional compact smooth 
manifolds without boundary, then diffeomorphic to the unit circle. Let a and b be the bilinear forms 
given by: 

a(u,v) = / AuAvdx, (4.27) 
Jn 

b(u,v) = / y~] d a/ 3ud af 3V dx. (4.28) 
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We consider the following variational problem (which has possibly only a formal sense) 



Find u e G V such that, \/v G V 

a(u £ ,v) +s 2 b{u e ,v) = (f,v), 

where the space V is the "energy space" with the essential boundary conditions on Tq 



dn\r 



0}, 



(4.29) 



(4.30) 



where n, t denotes the normal and tangent unit vectors to the boundary T with the convention that the 
normal vector n is inwards Q. It is easily checked that the bilinear form b is coercive on V. Moreover, 
we immediately obtain the following result. For all e > and for all / in V' , the variational problem 
( 14.291) is of Lax-Milgram type and it is a self-adjoint problem which has a coerciveness constant larger 
than ce 2 , with c > 0. 

The equation on ft associated with problem ( 14.291 ) is: 



(1 + £^)AV = / on ft, 



(4.31) 



as both forms a and b give the laplacian. As for the boundary conditions on Tq, they are "principal" 
i. e. they are included in the definition of V, (5.4). As for conditions on Ti, they are "natural", 
classically obtained from the integrated terms by parts. Those coming from the form b are somewhat 
complicated; we shall not write them, as the problem with e > is classical. For e = these 
conditions (coming from form a) are: Au = = 0, on Ti. 

As a matter of fact, the full limit boundary boundary value problem is: 



A 2 u° = f on Q 
u = = 0, on r 

An = on Ti 
_d_ 

dn 



(4.32) 



An = OonTi. 



Let us check that the boundary conditions on T\ (i. e; the two last lines of (5.6)) do not satisfy 
the SL condition for the elliptic operator A 2 . Indeed, proceeding as in sect. 2, by formal tangential 
Fourier transform 



+ d 2 ) 2 u = 0. (4.33) 

which yields that 

v = (Ae-\^ X2 + Cx 2 e-^ X2 ) (4.34) 

(as well as analogous terms with +|£| instead of — 1£|, which are not taken into account as exponen- 
tially growing inwards the domain). Here, according to SL theory, x 2 is the coordinate normal to the 
boundary, after taking locally tangent and normal axes, (which do not modify the equation A 2 ). The 
(tangential Fourier transform of the) boundary conditions on Ti are: 

(-£? + d 2 )u = (4.35) 

and 

d 2 (-f 1 + d 2 )u = 0. (4.36) 
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It is immediately seen that the previous solutions (14.341 ) with C = and any A / satisfy both 
conditions (note that its laplacian vanishes everywhere, then it vanishes as well as its normal derivative 
on the boundary). So, the SL condition is not satisfied on Y\. 

Before going on with our study, we note that the limit problem (14.321 ) implies an elliptic Cauchy 
problem for the auxiliary unknown 

v° = An . (4.37) 

Indeed, system (14.321 ) gives in particular: 

Av° = f on n 

v° = on Ti (4.38) 
= 0onIY 

which is precisely the Cauchy problem for the laplacian. 

As mentioned in section [3j this is a classical ill - posed problem, and the solution does not exist 
in general. Oppositely, uniqueness of the solution holds true (uniqueness theorem of Holmgren and 
analogous, see for instance ICoHi62ll ). 

4.2 The heuristic integral approach 

The aim of this section is the construction, in a heuristic way, of an approximate description of the 
solutions u £ of the model problem in the previous section for small values of e. 

From the general theory of singular perturbations of the form (14.291 ). we know that our assumption 

a(y, v) 1 / 2 defines a norm on V, (4.39) 

is crucial. Indeed, when it is not satisfied, the problem is said to be "non inhibited". In such a case, it 
has a kernel which contains non vanishing terms and then, it is easy to establish that the asymptotic 
behaviour of the solution u £ of (14.291 ) is described by a variational problem in this kernel. The previous 
fact is not surprising as soon as we consider the following minimization problem, which is equivalent 
to KM . 

Minimize in V, 

a(u £ , u £ ) + e 2 b{u £ , u £ ) - 2(f, u 



2u,.e ^ n,*..*S ^ 



Indeed, when e goes to zero, the natural trend consists in avoiding the a-energy which occurs with the 
factor 1 and leaving the 6-energy which has a factor e 2 . 

Clearly, this is not possible when (14.391) is satisfied since the kernel reduces to the zero function. 
Nevertheless, in our case, a(v, v ) = implies Av = and, as v E V, the traces of v and |^ vanish 
on To, so that (14.391 ) follows from the uniqueness theorem for the Cauchy problem. This uniqueness 
is classical, but the solution u is unstable in the sense that there can be "large u" in the V norm (or 
in other spaces) for "small /" in the V norm (or in other spaces). It then appears that the same 
reasoning shows that for small values of e, the solution u £ will be precisely among elements with 
small a(u £ , u £ ), that is to say with small Au £ in I? . 

4.3 The T layer 

Let us now build such functions u £ G V with very small || Au £ \\ L 2. The main idea is to consider 
functions in a larger space than the space of functions v of V such that Av = (which only contains 
the function v = 0). The functions of this bigger space will not satisfy the two boundary conditions 
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on Tq that are satisfied by any function of V. Then we shall modify it in a narrow boundary layer 
along Tq in order to satisfy the two boundary conditions with small value of a-energy. 
More precisely, let us consider the vector space: 

G° = {v £ C°°(Q), Av = on n, v = on T }. (4.41) 

Remark 4.1. We observe that every function of G° satisfies one of the boundary conditions on Tq 
which are satisfied by any element of V. For simplicity, we have chosen v = on Tq, but we could 
choose the other one ^ = on Tq as well. On the other hand, the regularity assumption C°° is 
slightly arbitrary. Since, we will consider the completion of G° with respect to some norm, this point 
is irrelevant. 

Obviously, as the Dirichlet problem for the laplacian on Q, is well posed in C°°, the space G° is 
isomorphic with the space of traces on T\: 

{w E C 00 ^)} (4.42) 
the isomorphism is obtained by solving the Dirichlet problem: 

Aw = on Q, 

w = on T , (4.43) 
w = w on T\. 

In the sequel, we shall consider indifferently the functions w on f2 or their traces won Ti. 

In fact, the exact function u £ is a solution of (14 .31b . which we are searching to describe approxima- 
tively in order to define a space as small as possible (incorporating the main features of the solution) 
to solve the minimization problem. More precisely, according to our previous comments, we are in- 
terested in the "most singular parts" of u £ in the sense of the part corresponding to the high frequency 
Fourier components. As we shall see in the sequel, it turns out that these singular parts may be ob- 
tained by modification of the functions tiona boundary layer close to Tq; this layer is narrower when 
the considered Fourier components are of higher frequency; in fact, the layer only exists because we 
only consider high frequencies. This allows to make an approximation which consists in using locally 
curvilinear coordinates defined by the arc of To and the normal, and handling them as cartesian co- 
ordinates. Clearly, this approximation is exact only on the very To, but more and more precise as we 
approach of To, i. e. as the considered frequencies grow. 

Once the layer is constructed, we compute the a-energy of it, as well as the e 2 6-energy of the 
(modified) w function, in order to consider the variational problem ( 14.291 ) in the restricted space. 

Let us first exhibit the local structure of the Fourier transform of w close to Tq. According to 
our general considerations on the heuristic procedure, w may be considered (after multiplying by an 
appropriate cutoff function) of "small support" near a point Pq of Tq. Taking local tangent and normal 
cartesian coordinates y±,y2, we have, within our approximation, 

d 2 d 2 , 

w = OonR x (0,t), (4.44) 



for some t > 0. Taking the tangential Fourier transform, we obtain: 

F(wj)(€i,y2) = Ae l?l1 ^ 2 + fie'^ lly2 . (4.45) 

It is worthwhile defining the local structure of w in the vicinity of Tq using the "Cauchy" data w and 
&2W on Tq (note that the solution of the Cauchy problem is unique, so that the Cauchy data determine 
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the solution). As w vanishes on To, the local structure is then determined by d 2 w on IV Taking the 
tangential Fourier transform this gives: 



\ Tf d ™i \ sinh(|^i|y 2 ) ,„ 



We now proceed to the modification of w into w a in a narrow boundary layer of To in order 
to satisfy (always within our approximation) the equation coming from (14.311 ) for small e. Using 
considerations similar to those leading to (14.44b . this amounts to 

id 2 d 2 \ (2) 
hence the tangential Fourier transform reads 



dy\, 



\^\ 2 + ^Y~' H™ a ) = 0- (4-48) 



Consequently, T{w a ) should take the form 

F(™ a )(£l,y2) = (a + jy 2 )e^ y2 + ((3 + 5y 2 )e-^\yi . (4.49) 

The four unknown constants should be determined by imposing that w a and ^w" vanish for 
^2 = and the "matching condition" of the layer, i.e., out of the layer, we want Wj to match with the 
given function wj. Since >> 1, then \£i\y 2 » 1 means that y 2 » r^r (but we still impose 
that y 2 is small in order to be in a narrow layer of To); this is perfectly consistent, as we will only 
use the functions for large hence the terms with coefficients (5 and 5 are "boundary layer terms" 

going to zero out of the layer (i.e. for \y 2 \ » see perhaps BEc79ll or III191II for generalities 

on boundary layers and matching. This gives 



F[w 3 )(t l ,y 2 ) = f[-± ( -lfte-IW»). (4 . 5 0) 

This amounts to saying that the modification of the function Wj consists in adding to it the inverse 
Fourier transform of 

)(-V2e-^). (4.51) 



<%2 |y 2=0 y 

Defining on To the family (with parameter y 2 ) of pseudo-differential smoothing operators Sa(e, D\ , 
with symbol: 

Sa(e^ 1 ,y 2 ) = -y 2 e-^h(e^,y 2 ), (4.52) 

where h is an irrelevant cutoff function avoiding low frequencies; it is equal to 1 for high frequencies 
(see [EgMeSa07] for details), we see that the modification of the function w: 

5w = w a -w (4.53) 

is precisely the action of 5a(e, D\,y 2 ) on §^-(2/1, 0): 

5w = 5a(e,D 1 ,y 2 )—^(y 1 ,0). (4.54) 
oy 2 
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Let us now compute the leading terms of the a-energy of the modified function w a . 

Let v and w be two elements in G° and v a , w a the corresponding elements modified in the bound- 
ary layer. As the given v and w are harmonic in f2, the a-form is only concerned with the modification 
terms 5v and 5w. Then, within our approximation, we have: 

r r+oo 

a(v a ,w a )= dyx A{5v)A(5w) dy 2 . (4.55) 
Jt Jo 

To compute this expression, we first write v and w as sum of terms with "small support" (by 
multiplying by a partition of unity): v = T,jVj and w = TijWj. Then, within our approximation, 
the integral is on the halfplane R x (0, +00) of the variables y±,y 2 . Taking the tangential Fourier 
transform and using the Parceval-Plancherel theorem, we have 

a(v a ,w a ) = X j:k / d6 / ( -g)6a(e,Z,y 2 )F(-± ) x 
J -00 Jo v a V2 / \y 2 =0 

{^l-^ie^r^ )dy 2 . 
Hence, on account of (14.521) and integrating in y 2 , this yields 

a(v»,w a ) = Z j , k r°°2|ei|^ ^T* h 2 (e,^y 2 )d^. (4.56) 
J-00 oy 2 ly2=0 dy 2 ly2=0 

This expression (14.561) only depends on the traces _ (yi) and _ (yi)> which are func- 

tions defined on To- 

We now simplify this last expression using a sesquilinear form involving pseudo-differential op- 
erators. 

Indeed, denoting by P{^) the pseudo-differential operator with symbol 

p(£i) = (2|6l) 1/2 M^,y 2 ), (4.57) 

and summing over j and k, we obtain that 



^->=/ r /<l>il/<l>§f lr - 

4.4 Taking account of the perturbation term e 2 b. 

We now consider the minimization problem (14.401 ) on G° instead of on V. Obviously, the a-energy 
should be computed using formula (I4.58I ). This modified problem should involve the a-energy and 
the e 2 6-energy. A natural space for handling it should be the completion G of G° with the norm: 



Ml 



p{-)- 

ds dn\v 



2 



ds + b(v,v). (4.59) 

It is easily seen that G is the space of the harmonic functions of H 2 (Q) vanishing on To; according 
to (14.431 ) it may be identified with the space of traces H z / 2 {Ti). 

It will prove useful to write another (asymptotically equivalent for large |£i|) definition of this 
problem. Indeed, the elements w of G° (and then of G) may be identified (by solving the problem 
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(14.431) ) with their traces w on Ti. Moreover, as the functions w are harmonic, we may exhibit their 
local behavior in the vicinity of any point xq S Y\. Proceeding as in (14.441 ). (14.451 ) and taking only the 
decreasing exponential towards the domain (this is the classical approximation for the construction of 
a parametrix) we have: 

H*>)(£i> y*) = ^H(6)e- |€l|w , (4.60) 

where y\ , y% are the tangent and the normal (inwards the domain) vectors. Then, it is apparent that 
the 6-energy is concentrated in a layer close to Ti and we may compute it in an analogous way to the 
calculus that was done for the a-energy (14.58b - Indeed, using Parseval-Plancherel Theorem and within 
our approximation, we have 

+oo 



b(w,w) = dm ]T|<9 Q/ ^| 2 dy 2 (4.61) 

J-oo JO aj3 

hence, recalling (14.601 ) and integrating over yi, we get: 

/+oo 
\d\ 3 \^(w)\ 2 d^. (4.62) 
-oo 

Then, defining the pseudo-differential operator Q(-§^) of order 3/2 with principal symbol 

V2|£i| 3/2 , (4.63) 

or equivalently as previously: 

V2(l + |£i| 2 ) 3/4 , (4.64) 
we have (always within our approximation): 



b(v,w) = Q(—)vQ(—)wds. (4.65) 

We observe that the operator Q is only concerned with the trace on Ti, so that we may either write 
v, worv,w in (14.65b - 

The formal asymptotic problem becomes: 



Find v e £ G such that \/w G G 



f To P(^)P(§) ds + e 2 / ri Q{¥) Q(w) ds = (/, w) 



(4.66) 



4.5 The formal asymptotics and its sensitive behaviour 

In order to exhibit more clearly the unusual character of the problem, we shall now write (14.661) 
under another equivalent form involving only the traces on Ti. Coming back to (14.43b . let us define 
TZo as follows. For a given w S C°°(ri) we solve (14.431) and we take the trace of ^ on To, then 

tt~~ = n w. (4.67) 

on |r 
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Using the regularity properties of the solution of (14.431 ). it follows that 7£o^ is in C°°(ro). In fact, 
Rq is a smoothing operator, sending any distribution into a C°° function. Then, (14.661 ) may be written 
as a problem for the traces on Ti: 



Find v e e # 3 / 2 (ri) such that G # 3 / 2 (Ti) 

Jr P&ov £ P&^ ds + e 2 J Q(-§- s ) V c Q(JL) W ds = J a Fw dx, 



(4.68) 



where the configuration space is obviously H 3 / 2 (Ti). The left hand side with e > is continuous 
and coercive. We then define the new operators 

A = TZqP*P1Zo € C(H s (Ti),H r (To)),\/s,r G R, (4.69) 
B = Q*Q e C(H 3/2 (Ti),H^ 3/2 (T l ) (4.70) 

where TZq is the adjoint of TZq (which is also smoothing)), (14.681) becomes 

'A + e 2 B) v e = F, in H~ 3 / 2 (Ti ) . (4.71) 



Obviously, B is an elliptic pseudo-differential operator of order 3, whereas A is a smoothing (non 
local) operator. 

This problem is somewhat simpler than the initial one (as on a manifold of dimension 1), show- 
ing the interest of the formal asymptotics. It enters in a class of sensitive problems addressed in 
|EgMeSa07 ] section 2. It is apparent that the limit problem (for e = 0) has no solution in the distri- 



bution space for any F not contained in C°°. Indeed, on the compact manifold To, any distribution is 
in some H~ m (To) space, which is send into C°° by the smoothing operator A. 

Remark 4.2. The drastically non local character of the smoothing operator A follows from the fact 
that it involves IZq and TZq ( see A4.67i ). This is the reason why the problem may be reduced to another 
one on the traces on T\. The possibility of that reduction is a consequence of our approximation, 
where the configuration space is formed by harmonic functions. 
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